Quasi-static approximation of the interspike interval distribution of neurons driven by 

time-dependent inputs 
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Variability in neural responses is an ubiquitous phenomenon in neurons, usually modeled with 
stochastic differential equations. In particular, stochastic integrate-and-fire models are widely used 
to simplify theoretical studies. The statistical properties of the generated spikes depend on the 
stimulating input current. Given that real sensory neurons are driven by time-dependent signals, 
here we study how the inter-spike interval distribution of integrate-and-fire neurons depends on 
the evolution of the stimulus, in a quasi-static limit. We obtain a closed-form expression for this 
distribution, and we compare it to the one obtained with numerical simulations for several time- 
dependent currents. For slow inputs, the quasi-static distribution provides a very good description 
of the data. The results obtained for the integrate-and-fire model can be extended to other non- 
autonomous stochastic systems where the first passage time problem has an explicit solution. 
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I. INTRODUCTION 



The sequence of inter-spike intervals (ISIs) generated 
by a neuron is typically variable. This variability has 
two different origins. On the one hand, under realistic 
physiological conditions, the total input current driving 
each neuron varies in time, thus producing also a variable 
response. On the other hand, neurons also give variable 
responses when repeatedly stimulated with one identical 
signal [l|, [2j . This extra variability is grounded on several 
inherently stochastic biophysical processes in the intrin- 
sic dynamics of neurons and in synaptic transmission. In 
this paper, we study how these two combined factors (ex- 
ternal variable currents and internal random processes) 
determine the ISI distribution of a simple model neuron. 

We work in the limit of slowly varying signals. Many 
cortical neurons are driven by slow time-dependent stim- 
uli. Pyramidal neurons often receive synchronous coher- 
ent input, easily measured in functional magnetic reso- 
nance imaging, and local field potentials (LFP). These 
input signals function as an external clock modulating 
the response of the cell. In some cases, the coherent in- 
put is oscillatory Q. Examples can be found in hip- 
pocampal place cells that receive massive theta input 
from the septum Q or in olfactory neurons modulated 
by the respiration cycle [f|. In visual cortical areas, the 
time-dependent signal modulating single neurons is more 
irregular, as shown by the LFP of primary sensory ar- 
eas jgj. Other examples of slow signals driving single 
neurons are given by the modulatory effects of certain 
neurotansmittcrs as dopamine, serotonin, noradrenaline 
and acetylcholine. These neuromodulators control the 
general level of arousal @, H[ ■ They affect the intrinsic 
neuronal properties of their target neurons in time scales 
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of seconds ji.g., much slower than those governing spiking 
dynamics [9( . Neuromodulation is particularly important 
in setting the rhythmic properties of small assembles of 
neurons controlling central pattern generators [l(J fill ]. 
As a final example, slow adaptation processes can also 
be described in terms of intrinsic transient currents with 
long time constants [H, , that either decrease (adap- 
tation) [3, [l5| or increase (reverse adaptation) [l6| the 
excitability of the neuron. The variable nature of these 
biologically relevant processes motivates the exploration 
of the response of a stochastic neuron driven by a slow 
time-dependent stimulus. 

We focus on the Perfect Integrate-and-Fire (PIF) 
model neuron. This is certainly a simplified model, con- 
taining no information of the ionic processes underlying 
spike generation. However, when working with dynam- 
ical neuron models, it is often convenient to waive the 
description of detailed biophysical phenomena, for the 
sake of simplicity. One such simplified approach is the 
Integrate-and-Fire (IF) neuron model (see [13, EU and 
references therein), which can be derived from more re- 
alistic approaches via the spike-response model [lj|. In 
an IF cell, the passive membrane integrator properties 
are described, whereas the whole spike-generation mech- 
anism (including afterhyperpolarization and resetting) is 
replaced by a discontinuous ad-hoc dynamical rule, ap- 
plied whenever the membrane potential reaches a fixed 
threshold value. This simplification is justified by the fact 
that the action potentials generated by a neuron are, in 
most physiological conditions, indistinguishable from one 
another. The neural code is based on the times of firing, 
and not on the shape of the spikes themselves. Thus, 
it makes sense to renounce to the description of action 
potential generation, and reduce spike firing to a point 
process. The PIF neuron is the only model neuron for 
which the exact ISI density function is known, for sta- 
tionary inputs. Our aim is to obtain its ISI distribution 
also for variable stimuli. 

In this paper we derive the shape of the ISI distribution 
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for the case of time-dependent driving, in the limit of slow 
stimuli. We show examples for linear, exponential, sinu- 
soidal and random Gaussian input currents, and compare 
the results of numerical simulations with the theoretical 
expression. 



II. QUASI-STATIC APPROXIMATION 

In this section we derive a closed-form expression for 
the ISI distribution of a PIF neuron, driven with an arbi- 
trary (integrable) slow input current and White Gaussian 
Noise (WGN). 



The noise £ (t) has zero mean, squared intensity 2D and 
no temporal correlations: (£(*)£(*')) = 2D8{t-t'). The 
continuous dynamical rule of Eq. ([I]) is complemented by 
a discontinuous process describing spike generation: once 
the membrane potential v reaches the threshold Vt, the 
neuron fires a spike and its voltage is reset to the value 
v r . The model is invariant under voltage displacements, 
so we can set v r = with no loss of generality. 

This particular stochastic model has an explicit solu- 
tion for the ISI distribution /(r|/i) of the interval r, when 
the neuron is driven with a constant current /i [2(J, HO] , 
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A. The starting model 



The first two moments of this distribution are 



Variability in neural activity in response to constant 
stimuli has been studied both theoretically 0, [2l| and 
experimentally t 22j. This variability has several origins. 
First, there is input variability. Although the experimen- 
talist may drive a given cell with a controlled stimulus, 
neurons are also affected by other inputs. In the cor- 
tex, each neuron receives synapses of tens or hundreds 
of thousands of other neurons. In sensory areas, to put 
an example, a large fraction of these afferents represent 
feedback connections from higher processing stages [23| . 
Their activity does not represent the external stimulus, 
but rather, a top-down regulation of the perception pro- 
cess. Their variability is presumably due to trial-to-trial 
fluctuations in the level of attention of the subject, in- 
terference with other processing modules, and memory 
or habituation effects. Although all these factors re- 
flect presumably relevant mental processes (and not just 
noise) , from the operational point of view, the resulting 
variability is outside the control of the experimentalist, 
and therefore is often modeled just as a random input 
current. A second source of variability is given by the 
stochastic nature of ion channels 0, [251 l2r| , synaptic 
release [13, HH or sensory transduction [23|. These are in- 
herently random processes, that govern neural dynamics 
at the molecular level. Their intrinsic unpredictability, 
based on thermic or quantal fluctuations, is also mod- 
eled as a random input term. 

We represent all sources of variability (irrespective of 
their origin) as an additive stochastic input current, mod- 
eled as White Gaussian Noise (WGN). In this context, 
we study a PIF neuron described by a voltage variable 
v, a capacitance C, driven by a positive external current 
/i(f) > and an additive WGN £(t). In the first place, 
we consider a constant input current fi(t) = ji = enst. 
The dynamical equations read 
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We can additionally rescale the model by setting 
v/v t — > v, [ijCv t — > /i, and D/C 2 vf — > D. In the 
new scale, the above equations still describe the sys- 
tem with vt = C = 1. In this new framework, voltage 
(v) is adimensional, and current (fi) has units of inverse 
time [T] -1 . Since there is no intrinsic timcscalc, the time 
unit is arbitrary. For definiteness, hereafter we consider 
[T] = ms. 



B. The case of a two-step current 

In the previous subsection we introduced the stochastic 
model, where the PIF neuron was driven by a constant, 
positive current and WGN. However, many interesting 
cases are not captured by this simple picture. Three dif- 
ferent situations arise in more realistic conditions: the 
deterministic part of the input current /i can vary as 
time goes by, the correlation function of the noise £(£) 
can dependend ont — t' (colored noise, with perhaps also 
a non-stationary intensity D(t)), or both. In this paper, 
we focuss on the first problem in the limit of slow inputs, 
and we mention how the approach should be extended, 
to also encompass the case of non-stationary noise inten- 
sity, D(t) [3l|. We therefore consider an arbitrary signal 
n(t) feeding into Eq. |T]). 

The simplest case consists of a two-step current (see 
Fig. [T]a). Our quasi-static approximation assumes that 
Eq. (O is valid inside each interval. The approximation 
holds when either the change in the current is small, or 
when each interval is long enough. 

In Fig. [T]a we show an example of a two-step current 
and a typical realization of the system of Eq. |T]). The 
interval with a larger input produces a stronger response 
(higher spike count). As shown in Fig. [TJb the two ISI 
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FIG. 1: (Color online) Two-step current, (a) Definition of 
the stimulation parameters (left) and a typical realization of 
the dynamics (right), (b) The histogram resulting from the 
simulation of the dynamical system of Eq. fljl is represented 
by the gray area framed by the stairs in black line (-Aerials — 
25000 and bin width = 0.2 ms), whereas the analytical result 
of Eq. (0 is shown in continuous line (red online) . Dashed and 
dot-dashed lines represent the stationary densities on each 
subinterval, /(t|/xi) and /(r|//2), respectively. Parameters: 
m = 0.1 ms" 1 , fi 2 = 0.25 ms" 1 , 71 = 150 ms, T 2 = 100 ms, 
D = 0.005 ms' 1 . 



density functions for the constant currents of each subin- 
terval (dashed and dotted lines) contribute to the total 
ISI distribution (gray shaded area framed by stairs) con- 
structed from the collection of all ISIs produced inside 
the entire stimulation period T. 

To estimate the ISI distribution in the complete stim- 
ulation period we assume that a given ISI of length r was 
either observed during the subinterval T\ or the subin- 
terval T 2 (see Fig. [TJa). This assumption constitutes an 
approximation, since in principle, an ISI can have its 
first spike in T\ and the second one in T 2 . If T\ and T 2 
are sufficiently long, however, then the single ISI falling 
astride the two intervals can be neglected, compared to 
the numerous ISIs that lie entirely in one step or the 
other. The quasi-static approximation consists in assum- 
ing that the conditional probabilities of finding the ISI r 
in each subinterval are equal to the stationary probabil- 
ities given by Eq. ©, with the corresponding values /i 
within each period (see Fig. [TJa). Notice that, again, 
this is not strictly true. Given that T\ and T 2 are finite, 
boundary effects in the first and last ISIs could in prin- 
ciple modify the distribution. 

In this context, we write 



f(r) = f(r\ f x 1 )-f(T 1 )+f(r\ f i 2 )-f(T 2 ). 



(4) 



In this equation, /(Tj) weighs the probability of draw- 
ing ISIs from the interval Ti. In addition, the probabili- 
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where (Ni) is the expected number of ISIs in the i-th 
subinterval, (Ni) = Ti/(n). Using Eq. © with v t = 1 to 
calculate (t<) we get 
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Replacing Eq. © in Eq. (f?|) we obtain the ISI distri- 
bution for the complete stimulation period, 



f(r) 
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Eq. ||7J) implies that the overall density is a weighted 
average of the densities in each subinterval. In Fig. [TJb 
we show the histogram obtained from the simulation of 
Eq. |T]) with the step-like input of Fig. [TJa, as well as 
the one obtained from the analytical expression, Eq. |(7J| . 
Clearly, both densities agree. 



Generalization to arbitrary time-dependent 
currents 



Based on the quasi-static approximation, the result 
of the previous subsection can be extended to an arbi- 
trary time-dependent input current. Suppose we make 
a partition of the stimulation period, {0, ti, tj, ijv}, 
as shown in Fig. f2J The smooth positive function /j,(t) 
can be approximated by a sequence of N constant strips, 
equal to the mean value of the function in each subinter- 
val, {/xi, ...,/ijv}- With this discrete partition, and 
proceeding as in the previous subsection, we obtain the 
ISI distribution valid for the complete stimulation period, 
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FIG. 2: (Color online) Generalization to variable currents. 
The time-dependent input is approximated by a step-like sig- 
nal. The time axis is divided into N intervals. In each interval, 
the original current is replaced by a constant value. 



where Atj = tj — tj-%. The partition can be inhomoge- 
neous Atj ^ Atj, for i j, and the functions f(r\fij) 
are the stationary probability density of Eq. (2}, for the 
particular values of the current within each subinterval. 
The continuum case is reached by making Atj — -> 0, 



each ISI. In this case, the input current since the last 
spike is essentially constant, and the approximation is 
expected to have a negligible effect. 

Equation (fTT)|) is valid for PIF neurons with additive 
WGN of constant squared noise intensity 2D. A similar 
argument can be used if also the noise intensity varies 
slowly in time. Since the noise intensity does not partic- 
ipate in the weighting factors, the extension of the pre- 
vious result to this more general case is straightforward. 
In this case, we obtain 



f(r) 
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dt 



In summary, we have derived the ISI distribution of 
a PIF neuron valid in an extended stimulation period, 
when the input current (both signal and noise) vary 
slowly in time. 



III. COMPARISON TO NUMERICAL DATA 



f(r) = 



Jo dt Jo 
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Replacing f[r\fi(t)} in Eq. (HJ) by Eq. ([2]), the quasi- 
static ISI distribution in the entire stimulation period is 



V4ttDt 3 / n{t) dt 
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dt. 



In the case of a two-step current, we averaged the dis- 
tributions /(r|/Xi) at each side of the jump, neglecting 
the single ISI that fell astride the two steps. In the case 
of continuous currents, in the limit of infinitesimal par- 
titions, all ISIs fall on two or more partitions. The ap- 
proximation, therefore, involves all ISIs. The procedure 
is justified when the stimulus varies only slightly inside 



We next test the validity of Eq. (|10|) for different in- 
put currents: linear, exponential, sinusoidal and Gaus- 
sian. We compare the analytical result with numerical 
simulations, to assess the validity of the quasi-static ap- 
proximation. We show that the numerical and analytical 
distributions agree, when the stimulus varies slowly. As 
the time scale of stimulus variations becomes comparable 
to the mean ISI, the two results begin to differ. 



Linear input current 



A linear input current is defined by fi(t) = A% + 
{A2-AX) ^ Parameters Ai and A2 control the strength 
of the input, whereas the variation rate is determined by 
the length of the stimulation period T. For fixed Ai, if 
T decreases, the current varies faster. 

The analytical expression given by Eq. (fTTjj) is inte- 
grable for linear currents, and the exact solution reads 
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Notice that the expression in Eq. (|12|) does not depend within the quasi-static approximation the ISI distribu- 
on the length of the stimulation period T. Therefore, 
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FIG. 3: (Color online) (a) Linear current case. Parameters of 
the input: Ai = 0.25 ms -1 , A 2 = 0.50 ms" 1 , T = 1000 ms, 
and D — 0.00125 ms -1 . Histogram (shaded gray area framed 
by the stairs in black line) constructed by simulating the dy- 
namical Eqs. (fTJ) with Nisi = 10 6 and bin width = 0.05 ms. 
Continuous line (red online): analytical result of Eq. (I12|l . 
Continuous line with circles (cyan online): distribution for a 
constant current Eq. (T2J), with fj, = (Ai + A-i)j1. (b) Ex- 
ponential input current. Parameters of the input: A\ = 
0.25 ms" 1 , A 2 = 0.25 ms" 1 , r e = 100 ms, T = 1000 ms, 
and D — 0.00125 ms" 1 . As in (a), the shaded gray area 
framed by the stairs in black line corresponds to the histogram 
obtained from simulations of Eq. (fTJ with JVisi = 10 6 and 
bin width = 0.05 ms. Continuous line (red online): theoreti- 
cal ISI distribution of Eq. (j 1 Op . Continuous line with circles 
(cyan online): the distribution for the asymptotic constant 
current /i = A\. 



tion only depends on the initial and final current values 
Ai and Ai. The total time T taken to stimulate the 
neuron is irrelevant. 

In Fig. [3]a we show the histogram obtained from the 
simulation of Eq. JT]), and we compare it to the probabil- 
ity density predicted by Eq. (|12p for a large stimulation 
period, T = 1000 ms. The initial value of the membrane 
potential in different trials was randomly chosen between 
and 1, and ISIs were collected between the first and 
the last spikes within T. The predicted ISI distribution 
agrees well with the real density. For comparison, we 
also show the density for a constant current equal to the 
arithmetic mean of A\ and Ai. Clearly, this distribution 
differs from the one obtained with a linear input current 
of the same mean. 



For small T, the input varies faster and the quasi-static 
approximation begins to fail. This occurs because of two 
reasons. In the first place, if at a given time t the typ- 
ical variation scale in the stimulus is smaller than the 
mean ISI, it is no longer true that the probability f(r\fj,) 
is well approximated by a constant current. In the sec- 
ond place, as the time before the first spike and the time 
after the last spike become comparable to the total inter- 
val T, boundary effects are no longer negligible. In order 
to circumvent the impact of boundaries, the analysis of 
the failure of the quasi-static approximation is delayed 
to the case of sinusoidal stimulation (see below) where 
boundary effects are negligible. 



2. Exponential input current 

The exponential current is defined by fj,(t) = A\ + Ai ■ 
exp (— t/r e ). Again, parameters A\ and Ai control the 
strength of the input, while r e defines the characteris- 
tic decay time. In this case, the explicit solution of the 
quasi-static expression Eq. (jlpp is not possible, so we re- 
sort to its numerical integration (Simpson rule). 

In Fig. [3jb we show the results for a large character- 
istic time r e = 100 ms. The total stimulation period T 
was selected to be ten times the decay time r e . For even 
larger T, the ISI distribution is dominated by the asymp- 
totic value of the input. Again, the distribution obtained 
from the numerical simulation of Eq. {1} is in excellent 
agreement with the theoretical result of Eq. (|10|) . We also 
show the density for the asymptotic constant current Ai . 
The mismatch between this density and the real distri- 
bution shows that even for T — 10r e , the whole stimulus 
history influences fir). 



3. Sinusoidal input current 

The sinusoidal current is defined by fi(t) = A\ + Ai ■ 
sm(ujt). Numerical integration of Eq. (fT0|) was used to 
obtain the quasi-static approximation for the ISI density. 
In this case, the characteristic time of variations in the 
stimulus is given by w -1 . Due to the periodicity of the 
input, the total stimulation time can be made arbitrarily 
long, thus minimizing boundary effects on the ISI distri- 
bution. Therefore, in this case, we can unambiguously 
explore the validity of the quasi-static approximation as 
the stimulus varies faster. 

In Fig. 31a we show the quasi-static ISI distribution 
and we compare it to the histogram constructed from 
simulations for a slow driving frequency / = 10 Hz 
(lo = 207T s" 1 ), and a large total stimulation period equal 
to a finite number of cycles (T = 1000 ms). Clearly, both 
distributions agree. For comparison, we also show the 
density for a constant input, with /j, = A\. Notably, the 
density for the mean value is completely different from 
the real distribution. 
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FIG. 4: (Color online) Sinusoidal input current. Param- 
eters of the input: Ai = 0.50 ms -1 , A2 = 0.10 ms -1 , 
T = 1000 ms, and D = 0.00125 ms" 1 . Histograms (shaded 
gray area framed by the stairs in black line) are constructed 
by simulating the dynamical Eqs. {TJ, with ATjbi = 10 6 and 
bin width = 0.025 ms. Continuous line (red online): theo- 
retical ISI distribution of Eq. l|10f> . Continuous line with cir- 
cles (cyan online): distribution corresponding to a constant 
current of intensity fj, = Ai. (a) Low frequency stimulus, 
u — 207T s" 1 (frequency = 10 Hz), (b) High frequency stimu- 
lus, lo = IOOOtt s _1 (frequency = 500 Hz). 



FIG. 5: (Color online) Gaussian input current, with cut- 
off frequency f c . Parameters of the input (temporal do- 
main): Mean = Ai = 0.50 ms -1 , Standard Deviation = 
A 2 = 0.10 ms -1 , Total stimulation period = T = 1000 ms. 
(a,b) Two signals with different cutoff frequencies were used: 
f c = 50 Hz and f c = 500 Hz (see insets at top of each panel) . 
Histograms (shaded gray area framed by the stairs in black 
line) are constructed by simulating the dynamical Eqs. {l}, 
with JVisi = 10 6 and bin width = 0.025 ms. Continuous line 
(red online): theoretical ISI distribution of Eq. (|10[) . 



As the input frequency increases, the quasi-static ap- 
proximation Eq. (fTO)) begins to depart from the real ISI 
distribution. In Fig. 0Jb we show a faster sinusoidal cur- 
rent with / = 500 Hz (w = IOOOtt s" 1 ). With this input, 
the real and the quasi-static densities are noticeable dif- 
ferent, and actually the real distribution gets closer to the 
mean constant input density. Crudely, a fast oscillatory 
current produces a similar effect to that of a constant 
mean current, with higher noise amplitude. 



4- Noisy Gaussian input currents 



In particular, we tested two stimuli with different cutoff 
frequencies: 50 Hz and 500 Hz, respectively. As shown 
in the insets of Fig. [5j the signal with the higher cutoff 
frequency displays more rapid fluctuations. In Fig. [5ja 
we see the ISI distributions from both the simulation of 
Eq. ((T|) and the quasi-static approximation Eq. (flTJ|) , for 
f c = 50 Hz. The agreement between the two distributions 
is excellent. As the signal fluctuates more rapidly, the 
quasi-static approximation for the ISI distribution starts 
to deviate from the real behavior. Fig. Ob shows the 
distributions for a stimulus of f c = 500 Hz. In this case, 
the quasi-static approximation is no longer applicable. 



Finally, we test the quasi-static approximation with a 
more irregular input, as done previously for example in 
f32| . In our case, the stimulus consists of a Gaussian 
signal with a sharp cutoff frequency f c . Similarly to the 
case of sinusoidal current, we denote the mean and the 
standard deviation of the signal in the temporal domain 
by A\ and A2, respectively. These parameters control 
the strength of the input. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper, we use a quasi-static approximation to 
derive the ISI distribution of a PIF neuron driven by 
a time-dependent input current. Different temporal de- 
pendencies were used to test the approximation. Linear 
currents were employed first, since they allow an ana- 
lytic integration of the ISI distribution. Currents with 
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other temporal dependencies required the numerical in- 
tegration of Eq. (jlOp . Exponential input currents are a 
realistic description of adaptation currents in the brain 
(l3| . Sinusoidal stimuli represent periodic collective os- 
cillations, as for example, the theta rhythm in hippocam- 
pus. Random low-passed Gaussian driving is an approx- 
imation to the LFP in sensory cortices. In all cases, for 
slowly varying stimuli, an excellent agreement between 
the theoretical ISI distribution and the one obtained by 
simulating the dynamical equations was achieved. As 
expected, when the stimulus varies in time-scales com- 
parable to the ISIs themselves, the approximation begins 
to fail. 

In previous studies, sinusoidal input currents were clas- 
sified as either endogenous or exogenous [33| . Exogenous 
stimuli are external signals, that do not depend on the ac- 
tivity of the neuron under study. Endogenous inputs, in 
contrast, are instantaneously reset to a particular phase, 
every time the cell generates a spike. Endogenous sig- 
nals arc more convenient for mathematical manipulations 
[33 |. However, as discussed in (35[, they do not consti- 
tute a realistic description of external sensory or synaptic 
input. They can only be justified in the subthreshold, 
noise-driven regime, where spike generation is only pos- 
sible in specific stimulus phases. In this regime, neuronal 
dynamics is characterized by a typical time scale, associ- 
ated to the average time required by the input noise to 
drive the cell across threshold. When this time scale is 
an integer multiple of the period of the sinusoidal input 
current, stochastic resonance and multi peaked ISI distri- 
butions are obtained [H, M, M, M, H, M, HI 5H, El ■ In 
this paper, however, we have focused on positive currents 
fi(t) > 0, which mimic the supra-threshold regime of 
leaky integrate-and-fire neurons. In this regime, the drift 
current term /i(t) is sufficient to generate firing. Input 
noise randomizes the process of spike generation, intro- 
ducing spike-time jitter, but does not determine a char- 
acteristic waiting interval. Therefore, there is no noise- 
dependent escape time and the multi peaked ISI distri- 
butions characteristic of stochastic resonance are not ob- 



served. 

The basic formulas of this paper Eqs. (jTTJJ) to (|12[) are 
only valid for the PIF neuron. However, extensions to 
other neuron models can be developed. To that end, one 
needs to know how (r^) depends on the input strength fi 
(and possibly D), to be able to derive a relation analo- 
gous to Eq. ©. In addition, /(r|/ii) must be known, in 
order to further proceed to Eq. (fTtl]) or Eq. (fTTj) . There- 
fore, the approach we pose in this paper could also be 
applied in other areas. For example, in an impressive 
modelling work about the nature of the ion channel gat- 
ing, Goychuk and Hanggi developed a stochastic micro- 
scopic description of the transition from closed to open 
conformations of potassium channels [43| . In their model, 
the authors posed the problem in the Fokker-Planck for- 
malism and argued for appropriate boundary conditions. 
They explicitly solved the residence time distribution in 
the Laplace domain and obtained the "unusual" opening 
rate of this kind of channels. One of the key points in 
their analysis is that the activation rate depends on a 
potential (describing the conformational energy barrier), 
which in turn depends on the voltage. However, in real 
conditions this voltage depends on time, since the ex- 
tracellular environment surrounding any central neuron 
is modulated by the collective LFP. The Fokker-Planck 
formalism used in [43] can also be phrased within the 
Langevin approach used in our work. Thus, our method 
could be useful to extend their results to time-modulated 
potentials. This example shows the broad applicability 
of the quasi-static formulation. 

The knowledge of the ISI distribution obtained for 
time-dependent stimuli is useful to assess the effect of 
an evolving current in the firing dynamics of a neuron. 
This distribution characterizes the running rate code, as 
received by a downstream target: if the post-synaptic 
neuron has a fairly long integration time, the statistics of 
ISIs received during the integration interval depend on 
the evolution of the stimulus inside that period. In ad- 
dition, in the quasi-static limit, the ISI distribution can 
also be used to calculate mutual information rates [271 ] . 
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